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Abstract 

We study the region inside the event horizon of charged black holes in five dimen- 
sional asymptotically anti-de Sitter space, using as a probe two-sided correlators which 
are dominated by spacelike geodesies penetrating the horizon. The spacetimes we in- 
vestigate include the Reissner-Nordstrom black hole and perturbations thereof. The 
perturbed spacetimes can be found exactly, enabling us to perform a local scan of the 
region between the inner and outer horizons. Surprisingly, the two-sided correlators we 
calculate seem to be geometrically protected from the instability of the inner horizon. 



1 Introduction and Summary 



Generic black holes, those with non-zero charge and/or angular momentum, have a causal 
structure very different to that of the more familiar Schwarzschild black hole. This is a source 
of new and interesting phenomena. In addition to the outer event horizon, the geometry 
contains an inner horizon surrounding a singularity, which is timelike. This inner horizon 
is also a Cauchy horizon - it bounds a region outside of the Cauchy development of an 
initial data set. Finally, the inner horizon is a surface of infinite blue shift - analysis at the 
linearised level shows that an arbitrary perturbation gets magnified by the geometry near 
that surface [H 12131. 

It is then generally expected that the inner horizon is unstable to small perturbations, 
which back-react strongly on the geometry, likely resulting in a singularity. The spacetime 
effectively ends at or before one reaches the Cauchy horizon, and so the cosmic censorship 
conjecture [3] is upheld. However, the precise mechanism for this expected instability, and 
the resulting spacetime structure, have not yet been settled. For recent efforts towards 
clarifying this issue within classical general relativity see [Hj, and for a recent investigation 
of three dimensional rotating black holes in this context see [H] . 

This set of issues could be studied within the AdS/CFT correspondence [Zl [Hj although, 
as for other interesting phenomena 1 associated with black holes, to study the physics one 
would have to probe the geometry behind an event horizon. Early efforts to do that are fUl 
[TTJ [TJJ HB] • Here, we will use the set of ideas introduced in [H] and further developed in 
[THJ HEl HZl UH1 El EE EHJ- Namely, we will study a set of observables which depend in a 
precise sense on the physics behind the horizon 2 . 

The observables are correlation functions of operators which are inserted on the two 
disconnected boundaries of the eternal black hole spacetime [HJ. These are the objects we 
will use to probe the above questions associated with generic black holes. In the limit of 
large operator dimension, the correlators are approximated by the contribution from a single 
spacelike geodesic, which penetrates the event horizon. For a precise explanation of this 
approximation we refer the reader to [TH I17j. 

This is an explicit realisation of work by Israel 22\, in which parallels were drawn between 

the concepts of the thermofield formalism of field theory (see, e.g., |2HJ for a review) and black 

hole physics. The fictitious Hilbert space in the thermofield formalism is just the space of 

states behind the event horizon of the eternal black hole. By tracing over the fictitious states, 

we recover the thermal behaviour of the physical states. In the AdS version, correlators of 

1 See, however, recent efforts to construct naked singularities in asymptotically AdS spaces 0- 
2 Though, in the spirit of black hole complementarity should also be recoverable from observations 
outside the horizon only. 
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operators inserted on a single boundary will be ordinary thermal correlation functions, but 
one may also consider correlation functions of operators inserted on two opposite boundaries. 
As explained in , one can compute such correlators from those of operators inserted 

on a single boundary, by analytic continuation, simply shifting the argument of one operator 
by half a period of imaginary time. 

The outline of the paper, and a summary of results is as follows: 

We begin, in section 2, by introducing the five dimensional AdS-Reissner-Nordstrom 
(AdS-RN) geometry, following |24j . We discuss coordinate systems and complexification of 
the metric, and draw the Penrose diagram which encodes the causal structure. We should 
note that we will only be concerned with the non-extremal black hole, since the extremal 
solution has only one asymptotic region, so there are simply no such objects as two-sided 
correlators. This makes sense, since the extremal black hole is at zero temperature [21]. 

In sections 3 and 4, we investigate correlation functions of neutral and charged operators 
respectively. The discussion does not involve the subtleties of |17j . as there is no branch 
cut or singularity, due to the drastically different causal structure. We find, however, an 
interesting non-analytic feature of the correlation function, which exists only for black holes 
sufficiently far from extremality. It would be interesting to discover the meaning of such a 
phase transition in the field theory language. 

As emphasised in ^H] > it seems that the field theory only encodes the region of spacetime 
within the Cauchy horizon. In that sense, the AdS/CFT correspondence would uphold the 
cosmic censorship conjecture, at least for the observables we consider. It would be interesting 
to see if this is a more general feature of the field theory. 

In section 5 we turn to an analysis of a perturbed AdS-RN spacetime. For perturbations 
which move at the speed of light (ingoing or outgoing only), one can use the methods of 
Poisson and Israel [23] to find the perturbed metric exactly. By varying the time at which the 
metric is perturbed, one can perform a scan of (at least some part of) the region between the 
inner and outer horizons 3 , again using the two-sided correlators as a probe. We demonstrate 
this through the simplest example of an infinitely thin perturbation, an outgoing thin shell 
of null matter. 

We confine ourselves here to calculating the dependence of the correlators on the time 
at which the perturbation leaves the boundary. This demonstrates that the gauge theory is 
sensitive to a perturbation of spacetime which is localised purely behind the outer horizon. 
It would be interesting to read off from these correlators the dramatic behaviour expected 
when mass inflation [2H] sets in. We leave this to future work [26 , but we comment that 
so far it seems that our correlators are screened from this catastrophic instability by some 
3 This is somewhat similar to the discussion in |T2], though here the perturbed metric is on-shell. 
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subtle geometrical effects, which we explain below. 



2 Classical Geometry 

In this section we discuss the classical geometry of the electrically charged AdS-RN black 
hole, its complexification and causal structure. We restrict attention to aspects of the ge- 
ometry which are relevant for us, more details can be found in |24j . 

The Spacetime 

The five dimensional solution 4 takes the form [21] 

d. 2 = - /( r)dt 2 + *L + rMnl A = + $) dt (2.1) 

where 

f(r) = 1 + ^- ^ + ^ = ^, A(r) = r« + ZV - l 2 Mr 2 + l 2 Q 2 , (2.2) 

c = 2/v / 3, $ is a constant and M and Q are respectively proportional to the mass and 
charge of the black hole. The spacetime is asymptotically AdS, with curvature radius Z, the 
boundary theory living on K x S 3 . From now on, we will take 1 = 1. The boundary is at 
r = oo and there is a timelike singularity at r = 0. 

Horizons of the metric are given by the real positive roots of A(r). We will often work 
with the variable x = r 2 , A(x) being a cubic in x. Since, with M > 0, 

A(0) = Q 2 > and ^(0) = —M < 0, (2.3) 

dx 

one root of A must occur for x < 0. We denote this negative root by — Xq = — r 2 ,. There are 
then at most two positive real roots, denoted by x± = r 2 , with x + > X-. We will thus write 

A(x) = (x + Xq)(x — X-)(x — x + ). (2.4) 

It is easy to show that xq = 1 + x_ + x + , so we will only need to specify x±. The non-extreme 
black hole, with x + > x_, has an outer event horizon at r + and an inner Cauchy horizon at 
r_. The extreme solution, with has a single event horizon at r + . 

The surface gravities at the outer and inner horizons are, respectively, 

(x + x±)(x + -X-) . . 

k± = §75 > (2.5) 

x± 



4 A solution of five dimensional Einstein-Maxwell theory with a negative cosmological constant. The ten 
dimensional origins are discussed in |24j. 
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the temperature of the black hole being T = k + /(2tt) = 1/(3. The dual field theory is thus 
in a thermal state at temperature T. One can then work [21| with either fixed charge (the 
canonical ensemble) or fixed potential (the grand canonical ensemble). In the latter case, the 
electrostatic potential in the field theory is given by the constant $, which can be fixed as 
$ = (l/c)(Q/x + ) such that A t {r + ) = 0. In terms of the bulk physics, $ is the electrostatic 
potential difference between the horizon and infinity. 

Embedding into the Complex Plane 

There are six distinct regions of the spacetime, the Penrose diagram consisting of an infinite 
sequence of these basic blocks. As was emphasized in [THj . and as will be discussed below, the 
region behind the inner horizon does not seem to play any role in the field theory, therefore 
we can restrict attention to one copy of these basic blocks. 

To describe the global structure, one can pass to Kruskal-like coordinates to describe 
the maximal analytic extension of the space. Alternatively, as in [16*1 ITTj IT%1 |E] , one can 
use different Schwarzschild coordinate patches to describe the global extension. They are 
embedded into complex Schwarzschild time, with t — + He having a constant imaginary 
part in each patch. This is shown in figure 1, in which we define = in the right asymptotic 
region. Then, as in [T7], crossing the event horizon at r + shifts t E by 7t/(2k + ) = i/3/4 so, 
to move from a point on the right-hand boundary to the symmetric point on the left-hand 
boundary, we shift t — > —(t + i(3/2). Crossing the Cauchy horizon at r_ instead would shift 
tE by 7r/(2/t„) although this will not be relevant to the boundary field theory 5 . 

The Penrose Diagram 

The basic features of the Penrose diagram are shown by the behaviour of radial null geodesies. 
As in the uncharged case [Hj , the boundaries and the singularities cannot both be drawn as 
straight lines. For radial null geodesies, we work with the Lagrangian 

C = -/(r)i 2 + J£j = (2.6) 

where dots denote differentiation with respect to the affine parameter A. There is a conserved 
energy, E, associated with the Killing vector d/dt, in terms of which the geodesic equations 
become 

i = #y - 2 = E 2 . (2.7) 



The embedding into complex time can be repeated infinitely in both directions, by shifting ts appropri- 
ately, as one crosses each horizon. 
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Figure 1. The embedding into complex Schwarzschild time, showing the constant imaginary parts of t = 
th + He- 



For ingoing geodesies which start at the boundary r = oo at t = 0, the time t as a 
function of r is given by 

f°° dr' r°°dx 

~ Jr W) ~ Jr* Y {X + X ){X - X^){X - X + )' ^ ' ' 

Explicitly, for r < r_, 

t(r) = —, % r ftan^ 1 (— \ - ^] - — tanh" 1 (— 

{x + x_){x + x + ) V V r o/ 2; k_ 



«+ V r + 

where the pure imaginary terms arise from integrating over the two poles at x — x±. These 
are just the shifts in ts discussed above. Such a geodesic will reach the singularity r = at 
a finite time 

3/2 • / 1 1 

7T X n Z7T / 1 1 



tsins 2 (x + z-)(so + x + ) 2 U+ «- ' ' (2 ' 10) 



whose real part is positive. Conversely, a geodesic which hits the singularity at tt = 0, must 
leave the boundary at time 



tr 



7T 



x, 



3/2 



2 (x + xJ)(xq + x + ) 



<o, 



(2-11) 



from which we conclude that one cannot draw both the boundaries and the singularities as 
straight lines in the Penrose diagram. Choosing the boundaries to be straight lines, then 
the singularities must be bowed out, the situation demonstrated in figure 2. However, it will 
become clear that this behaviour has no influence on the boundary theory. In particular the 
critical value t c in the Schwarzschild case of [Tj\ plays no role in the discussion of correlation 
functions here. 



t = 




Figure 2. The Penrose diagram repeats itself in both directions. Radial null geodesies are dotted red, radial 
spacelike geodesies are solid blue, the inner and outer horizons are dashed and the singularities are bold. 
We have chosen to draw the boundaries as vertical lines, in which case the singularities are bowed out. The 
values of t all denote those of t^, the real part of complex time. 
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Kruskal Coordinates 

As in the asymptotically flat case [27], we need two Kruskal-like coordinates patches to 
cover the region < r < oo, one being valid through the outer horizon, the other being valid 
through the inner horizon. In the usual way, we first define the tortoise coordinate 

r . = £ J* + C = 4tan- Q + J_ h (iZi) - J_ ln (^) ( , 12) 

where Xq = x / ((xo + x + )(xo + and we have chosen the constant C so as to make r* 
real for r > r + . As we cross the horizons from this region, we pick up the same constant 
imaginary terms as in the expression (J2.9)) for t, giving an embedding of the spacetime into 
the complex coordinate plane. 

To cover the entire range of r, define the light cone coordinates u = t — r*, t> = t + r*. 
For r > r_, transform according to 

U + = -e~ K+u = T + - X+ V + = e +K+v = T + + X+ (2.13) 

whereas for r < r+, take 

U~ = -e +K - u = T+ - X + , T/~ = e^"" = T~ + X~. (2.14) 

The outer horizon is given by T +2 — X +2 = and the inner has T~ 2 — X~ 2 = 0. The metric 
becomes 

ds 2 = -L/(r) e T2K±r *(-(dT ± ) 2 + (dX ± ) 2 ), (2.15) 
where for r > r_, r(T + , X + ) is determined implicitly by 



_ _^+^ _ _ e ^K+r» _ _ e ^K + x tan ^r/roj 21 2 (2 16) 

\r + r + J \r — r_ 

and for r < r + , r(T~, X - ) is determined by 

T -2 _ x -2 _ _ e -2K-r, _ _ e -2 K _^tan- 1 (r/r ) f T + + r \ 1 + f T-—T_ \ ^ (2.17) 



r i — r / \ r_ + r 



However, the boundaries and singularities are covered by different coordinate patches, which 
makes comparison of them impossible. 

One can instead follow [2TJ and use the coordinates defined in ()2.13|) for the regions 
r > r + and r < r_. In the region r_ < r < r+, we instead define the light cone coordinates 
u = t + r*, f = — (t — r*) and use 

[/ = e K+u = T-X, 1/ = e K+,; = T + X. (2.18) 
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The metric can be written everywhere as 

d S 2 = \l^A e -^+r, (_ dT 2 + dx 2) ^ 1Q j 

which is analytic except at the inner horizon r = r_. The radial coordinate r is determined 
implicitly by 

T 2 - X 2 = zpe^+^otan-V/n)) ( \ r ~ r +\ \ ( r + r - A +/ (2.20) 

\ r + r + / \ \r — r_| / 

the =F sign corresponding to the regions r > r + and r < r_, and r_ < r < r + , respectively. 

Now we can compare the boundaries and singularities in the first set of coordinates. In 
the limits 

r -> oo =>• T 2 - X 2 -> - e ™ +a; ° < -1, (2.21) 
r -> T 2 -X 2 ^-l (2.22) 

so in a (T, X) spacetime diagram, the hyperboli representing the singularities will be further 
away from the origin than the hyperboli representing the boundaries. This is opposite to the 
Schwarzschild case ^7]. Following those arguments, one sees that on the resulting Penrose 
diagram if we draw the singularity r = as a vertical line, the boundary r = oo must be 
bowed out. Alternatively, we can use a conformal transformation to make the singularity 
bowed out while keeping the boundary vertical. The latter case is shown schematically in 
figure 2. (In this case, to bring the boundaries of the spacetime to a finite coordinate distance 
we can, for example, let U = e 7TK+x ' ^ 2 tanu and V = e 7TK+x '°^ 2 tanv.) 

3 Neutral Correlation Functions 

We now discuss features of correlation functions of operators which are electrically neutral. 
When we put insertions on both boundaries, the correlation functions are dominated (in the 
limit of high dimensional operators) by spacelike geodesies, as was explained in [Tlj. Such 
geodesies are only sensitive to the metric, and not to the background electric field. In the 
present section, we can restrict ourselves to symmetric geodesies (which reach the boundaries 
at the same value of ti = t ), since all others can be obtained by time translation. 

We start in subsection 3.1 by discussing the qualitative features, before presenting a more 
quantitative calculation of the correlation functions. As in the null case above, the geodesies 
are specified by the value of the energy E. In subsection 3.2, we turn to discussing in detail 
the relation between E and the boundary time, t(E), and in subsection 3.3 we compute the 
length of the geodesies, which will determine the correlators in our approximation. 
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3.1 Qualitative Features 
Radial Spacelike Geodesies 

For radial spacelike geodesies, and normalizing appropriately, the Lagrangian is 



C 



-f(r)i 2 + 



f(r) 



1. 



The geodesic equations are 



t 



E 

W) 



E 2 + f(r) = E 2 - V(r) 



(3.1) 



(3.2) 



the latter describing a particle of unit mass and of energy E 2 moving in an effective potential 
V(r) = —f(r). General properties of the geodesies 6 can thus be read off from figure 3, in 
which the effective potential is shown. 

V(r) 




Figure 3. The effective potential in the non-extreme case is shown in red, and the extreme case in blue. 

The extreme case is relatively uninteresting since there is no region between the horizon 
and the singularity which geodesies can probe. Geodesies with E 2 > fall all the way into 
the singularity; those with E 2 = reach the single event horizon before returning to the 
(same) boundary. 

The non-extreme case has the following qualitative features: 

• All geodesies with E 2 < V max return to the asymptotic region r = oo, so will connect 
the two boundaries. The endpoints lie on the two opposite boundaries, as can be seen 
by inspecting t(X). 

6 The uncharged, non-rotating geodesies discussed here resemble geodesies in the AdS-Schwarzschild ge- 
ometry which carry angular momentum |17| . 



9 



Geodesies with < E 2 < V max will penetrate the outer horizon, and reach some 
minimal turning point between the two horizons. They cannot penetrate the inner 
horizon, and will in fact accumulate at some finite distance from it (the radius for 
which the potential attains its maximum V ma _ x ), as E 2 — > V m 



max- 



• Geodesies with E = V max penetrate the furthest before returning to the boundary. 
These geodesies will turn out to have diverging boundary times. 

• Geodesies with E 2 > V max will fall into the singularity, where our approximation breaks 
down. As the boundary time diverges before we reach this range of E 2 , it seems that 
such trajectories do not contribute to any gauge theory process. 

Apparently, at least for the observables we are interested in, the region beyond the Cauchy 
horizon is not encoded in the gauge theory. Similar conclusions about three dimensional 
rotating black holes have been obtained in ^H]- Moreover, the geodesies do not even reach 
the Cauchy horizon itself, reaching at most some minimal distance 7 from it. We will see later 
that, even when perturbing the spacetime, the observables of interest are screened from any 
dramatic behaviour associated with the Cauchy horizon by this geometrical effect. 

Rotating Spacelike Geodesies 

Let us add an angular momentum L, conjugate to the azimuthal angle <fi on the boundary 
three-sphere 8 . It is the conserved quantity associated with the Killing vector d/d<p, the 
relevant geodesic equation being 

j>=^. (3.3) 

In this case, the effective potential is modified to V(r) = — /(r)(l — L 2 /r 2 ) which, in addition 
to r + and r_, has an additional root at r = \L\. A sketch of one possibility is given in figure 
4. 

In the vicinity of the singularity at r = 0, the potential is very different to the non- 
rotating case: for an arbitrarily small angular momentum, it turns around at small values of 
r. Consequently, for sufficiently large E one can come arbitrarily close to the singularity and 
still return to the asymptotic region r = oo; naively the corresponding two-sided correlators 
are sensitive to the region inside the inner horizon. 

However, a closer inspection reveals that such geodesies are irrelevant for the two-sided 
correlators. Once again, as E 2 — > Vm ax , the boundary time diverges. In addition, the 



7 This minimal distance can be large for far from extremal black holes, and goes to zero in the extremal 
limit. 

8 For more general rotation, one replaces L 2 by the relevant Casimir operator. 
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V 




r 



Figure 4. The effective potential for a rotating geodesic with x~ < L 2 < x+. Other choices of the relative 
magnitude of L 2 , X- and x+ simply switch the roots. If L 2 is equal to either a;_ or x + , then those two roots 
will coalesce. The local maximum V max is also shown. 

behaviour of t(X) shows that these geodesies connect two boundary points on the same side 
of the Penrose diagram. It is not clear what role, if any, such components play in the field 
theory dual, since inner horizon instabilities are likely to change the causal structure and 
eliminate these additional boundaries. 

The qualitative details of the rotating trajectories with E 2 < V mSuX depend only on the 
relative values of r + and L. For \L\ > r + the geodesies do not penetrate the outer horizon, 
making them uninteresting for our purposes. On the other hand, for \L\ < r + the situation 
is basically unchanged from the non-rotating case, the geodesies for which < E 2 < V max 
being relevant to us. 

3.2 Boundary Time 

In this subsection we describe the relation between (the real part of) the boundary time to 
and the energy E. We are interested in symmetric geodesies which start at a point to o n the 
boundary, and end at the point —(to + if3/2) on the other boundary. Taking te to denote 
the turning point for a trajectory of energy E 2 , we have 



E 2 + f(r E ) 



0. 



(3.4) 



As in [T7], we must have ii(r^) = and ^(r^) 



(3/4. Such geodesies satisfy 




(3.5) 
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where V(r) is the effective potential described above, for either rotating or non-rotating 
geodesies. Since the qualitative features are similar, we will only be concerned with the 
latter, for which V(r) = —f(r). 

The two features prominent in the analogous discussion of [T7], for the AdS-Schwarzschild 
metric, are the branch cut at to = (corresponding to E — 0) and the existence of a critical 
value to = —t c (corresponding to E — oo), where the correlators naively diverge. The latter 
feature was interpreted as a signal of the singularity encoded in gauge theory correlation 
functions, the former making that encoding a subtle one. In our discussion these features 
are no longer present. The special values of the energy are E = and E 2 = V max , and we 
discuss them both below. 

For general energies, one can compute ()3.5|) in terms of elliptic integrals 9 , which we do in 
the appendix. The result will depend both on the roots of the cubic A(x) and on the roots 
of the other cubic appearing in the denominator of ()3.5j) . 



A (a;) = E x + (x + x )(x — X-){x — x + ) = (x + x )(x — xJ){x — x + ) 



(3.6) 



where the roots x + > a;_ > > — x and x + = r%. 

From the appendix, and dropping the pure imaginary piece, we have 



—E 



^X+(X- + Xq) 



(X + + X )(X_ + X )(X + X_)(X + X + ) 



n((j),a ,k) 



x 3 (x 



x. 



(x + — xJ)(x^ — x^)(x + x_)(x + — x_) 



Q> 2 _ 



,y 3 ( , >-» 
\Xj _i_ ^ -Xi _|_ 



X. 



(x + — x+){x- — x + ){xq + x + ){x + — X-) 

F(4>,~k) 



n(0,ai,A:) 



x 3 



(X- 



where 



~2 



x + + x 



X- + Xq J \X + + Xq 



4- Xq)(x. 

X_ + Xq 



xJ)(x- 



— X-\ 



(3.7) 



a 2 . 



x + + x 

X- + Xq 



X. 



X. 



x± x_ 



a 2 . 



Xq 



Xq 



Xa 



X a_ Xj 




sin 




This is a fairly unilluminating expression, but more information can be seen if we plot 
t(E) numerically, for various choices of black hole parameters. We find two distinct types of 
behaviour, shown in figure 5: 

9 We use the notation of j^Hl for the elliptic integrals. 

10 This expression is actually valid for any geodesic, not just symmetric ones. We simply have to replace t 
with ib — itpj these values denoting, respectively, the values of ti, at the boundary and at the turning point. 
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• The left-hand plot is representative of the behaviour for x + ~ x_, where the black hole 
is near-extremal. The boundary time t is then a monotonically increasing function of 
energy E. 

• The right-hand plot, on the other hand, represents the behaviour for x + » x_, this 
limit being one of large mass and small charge. 

• In both cases, we clearly see a divergence in t as E 2 — > Vm ax . 




Figure 5. The left-hand plot shows t{E) for x+ ~ X- whereas the right-hand plot has x+ >> X-. The 
dotted lines indicate E 2 = Vmax- 

There is thus a range of parameter space in which, for each value of boundary time t, 
there exists a unique geodesic connecting the two boundary points (the left-hand plot in 
figure 5). However, outside of this range (the right-hand plot in figure 5), there are values of 
boundary time for which there are multiple geodesies connecting the two boundaries, each 
with a different value of E. Which one dominates the path integral will depend on the 
relative proper lengths. 

Since the cross-over between different saddle points will give rise to non-analyticity of 
the correlation function, we see that the existence or otherwise of non-analytic behaviour 
depends on the black hole parameters, x±, in an interesting way 11 . We comment further on 
this feature below (in subsection (3.3)), when discussing the proper length of these spacelike 
geodesies. 

11 Similar changes of behaviour as a function of parameters was noticed and interpreted in [22], and it 
would be interesting to relate these two types of phase transition. 
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The precise point of cross-over between these two types of behaviour, as a function of 
parameters, is difficult to pin down analytically but is possible to analyse numerically. For 
example, figure 6 shows a region of x± parameter space, the circles denoting that t is a 
monotonically increasing function of E for those choices of parameters, and the crosses 
denoting that t has a turning point 12 . 

x + 



10 + 



X. 



Figure 6. This plot shows how t(E) behaves as a function of the parameters x±. Circles denote points for 
which t is a monotonically increasing function of E, and crosses denote points for which t has a turning point 
for some value of E. 



Small E Geodesies 

It is easy to see that E = corresponds to to = 0, both for rotating and non-rotating 
geodesies. For E = 0, we have 

i = t(A) = 0, r 2 = /(r), (3.9) 

so these geodesies start at the boundary, pass straight through the Penrose diagram and end 
at the other boundary, without ever penetrating the horizon. 

To derive the small E behaviour of ([3.7)1 . we need the ^-dependence of the roots xq,X- 
and x + . This is determined by the following set of equations: 

Xq — X_ — X^. = Xq — X— — -|- E , 

X — Xj^ — XqX^ — XqX— = — XqX^. — XqX_, (3.10) 

XqX-X + = XqX-X + . 

12 It is curious that the boundary between these two types of behaviour is a straight line. 
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Since we are interested only in small E, we can solve these by taking 



x + (1 + a- + a+)E 2 + 0(E 4 ), 
x^ + a^E 2 + 0(E A ), 
x + + a + E 2 + 0{E A ) 



(3.11) 



where 




(3.12) 



It is important that the possible order E terms in (|3.11j) necessarily vanish, as long as 
x + ^ which in turn implies that all odd powers of E also vanish. We then have 



= + 1 E 2 + C(E 4 ), k = k + k 1 E 2 + 0{E i ), a 2 = a 2 + a\E 2 + 0{E 4 ) (3.13) 



(for each of the as), and where the quantities without a tilde are independent of E, being 
functions of the roots xq, X- and x + only. The elliptic integrals can then be Taylor expanded 
in each of their arguments using the well-known expressions for the various derivatives [28J. 
The net result is that 



for some non-vanishing constant t±. This expression for t(E) is similar to the finite mass 
black hole of [T7j . and in particular is an analytic function near E = 0. 

The Regime E 2 — > V max 

The integrand in ()3.5|) is divergent at the turning point of the trajectory, where E 2 = V(r). 
This singularity is integrable for E 2 ^ V max , leading to finite boundary time to for < 
E 2 < V max . On the other hand, when E 2 = Knax, the two positive roots of A(x) coincide 
(x+ — X-), leading to a logarithmic divergence in ()3.5|) . We emphasise that this is true for 
both the rotating and non-rotating geodesies. 

The divergence is shown clearly in figure 5, but we can see it analytically from the exact 
result ()3.7|) in the following way. Take e and note that k and each a 2 go to one 

faster than goes to ir/2, since 



to = hE + 0(E 3 ) 



(3.14) 



k = 1 + 0(e) 



a 2 = l + 0{e), 4> = - + 0{e 




(3.15) 



Then note that 




(3.16) 



X-L. — X — — 



and use [2H] 



n(0, 



i,k) 



(i - k 2 ) 



i 
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The problematic (1 — k 2 ) factors cancel, giving 

x 5/2 

to = -E— — _ . — In(tan0 + sec0) 

[X + + Xo) 6 ' [x+ — X-)[X + — x + ) 

~ -ln(cosvr/2) + finite, (3.18) 
since x_ < x + < x + . 

Since the boundary time diverges as E 2 — > V max , the gauge theory does not seem to 
encode the regime E 2 > V^g^. 

3.3 Correlation Functions 

The proper length of both the rotating and non-rotating geodesies is given by 

/*f max J ~, 

L = 2 - (3.19) 



^E 2 - V{i 



with V(r) being the relevant potential in each case, is the same turning point as in 
the previous subsection and, since we will only be concerned with the non-rotating case, 
it is again given by solving ()3.4j) . The upper limit, r max , is a long-distance radial cutoff, 
dual to the UV cutoff in the gauge theory. As we take r max — > oo, the integral diverges 
logarithmically. To regularise it we subtract the divergent piece arising in the pure AdS 
case (obtained by setting Xq — 1 and x_ = x + = in the above expression). This standard 
process is dual to renormalisation of the boundary theory [2*H] . 

E = Geodesies 

We start with the E = symmetric geodesies, which approximate the boundary correlators 
with to = 0. The proper length of such a geodesic is 



dr /""max / 7f 

L = 2 I -f== / dxj- . . . . (3.20) 

V/(r) Jr% y {x + x )(x - x_)(x - x + ) 



where the turning point for E = is r = r + , the location of the outer horizon. We can 
compute ()3.20|) in terms of elliptic integrals, giving 



a/x+(x_ + Xoj 

where 



[ar_ F(0, k) + (x+ - x_) U(<f>, x + k 2 /x^, k)] (3.21) 



sin- 1 J[ X - + XQ l^ ~ X j, k=J X -( X+ + Xo) =,[^±. (3.22) 



(x+ + x )(r max -x_)' ' \x + (x_+x ) y ft- r_' 
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In the limit r max — > oo, the logarithmic divergence appears in 11(0, x + k 2 /x-, k) (whereas 
F(<j),k) remains finite). To extract this divergence, we note that a 2 = x + k 2 /x- > 1, so we 
can use [30J 

.2 ,x . TT,, ,2,2 ^ , t,,a ,.\ , 1 i_ /A(0)+ ft tan< 



n(0, a 2 , fc) = -11(0, fc 2 /a 2 , fc) + F(0, fc) + — In - " ^ ^ (3.23) 

zpi \ A(0J — pi tan i 



where p\ = \J (a 2 — 1)(1 — k 2 /a 2 ) and A(0) = a/1 — k 2 sin 2 0. Expanding to first order in 



L = ~ = [x+ F((f),k) - (x+-x-) Il((f),x-/x + ,k))+\n ( — — j+ln^J, 

y/x+(x- +x ) \x + x„+x + J 

(3.24) 

the logarithmic divergence being manifest. The renormalised result is thus 

2 ( 4 

[x + F(0, fc) — (x + — X-) 11(0, X-/x + , k)] + In 



y/x+(x- +a; ) " \x + x_ + x +i 

(3.25) 

where now 

< E < V m£LX Geodesies 

As in [T^J, as we increase the energy, the geodesic will penetrate some distance inside the 
horizon. The proper length of the (non-rotating) geodesic is now: 

2 

rTxaax /* r max I ^ 

C ~ 2 Jr E y/E 2 + f(r)~ J r % E 2 x 2 + {x + xo) {x-x_){x- x + ) 1,2,1 

which will depend on the roots of the cubic A (a;), defined in (|3.6j) . Formally, the result is as 
in ([3.25)1 . but with the ^-dependent roots: 



yjx+{x- + x ) 



x + F(0, k) — (x + — xJ) 11(0, x_/x + , k) 



+ In 



(3.28) 
where 

= sin" 1 , l^±*°, k = , + (3.29) 

To take the small E limit of this expression, we proceed as in the previous subsection, 
using the expressions (|3.1ip . The upshot is that 

L(E) - L(E = 0) = L\E 2 + 0(E 3 ) (3.30) 
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for some non- vanishing constant L\. Using the previous expression (|3.14|) for the small E 
behaviour of the boundary time, we see that 



L(E) - L(E = 0) ~ t 1 



(3.31) 



The result is therefore analytic around E = 0, the branch cut structure found in for the 
infinitely massive black hole being absent here. It is more similar to the finite mass black 
hole (m. 

On the other hand, as for the computation of the boundary time above, the integral ()3.19|) 
will have an additional logarithmic divergence at the lower limit for E = Vm ax , which we 
can again extract analytically, or observe numerically. In figure 7, we plot L(E), the proper 
length as a function of energy, for the two different regimes of the previous subsection. The 
parameters are the same as those in figure 5, the left-hand plot being representative of the 
behaviour for x + ~ x_, and the right-hand plot being representative of the behaviour for 
x + » x_. The divergence as E 2 — > V mSQi can clearly be seen. 



L - L(0) 



L-L(O) 



E 




Figure 7. The left-hand plot shows L(E) — L(E = 0) for x+ ~ x_ whereas the right-hand plot has x + >> x- 
The dotted lines indicate E 2 = V ma ^- 



Non- Analytic Behaviour 

We have found that the subtle features of the infinitely massive AdS-Schwarzschild black hole 
are absent in our calculation. This is not surprising, since even very far from extremality we 
are expected to reproduce only the features of the finite mass black hole. 

On the other hand we find a new source of non-analytic behaviour which is much less 
subtle, coming from multiple real geodesies connecting the same boundary points. The 
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precise type of non-analyticity depends on the parameters of the black hole. The possibilities 
are shown in figure 5, the plots of t(E). 

For large enough boundary time, there is always a unique geodesic connecting the bound- 
ary points. As we decrease the value of the boundary time, there is a critical time t crit (and 
corresponding energy E crit ) for which additional geodesies appear. This is seen in the right- 
hand plot of figure 5, for the case of a black hole far enough from extremality. The type of 
non-analyticity then depends on which geodesic has the least action, which can be read from 
figure 7, the plots of L(E). To make this process more transparent, in figure 8 we plot L(t), 
the value of the correlator as a function of boundary time for all branches, and for different 
values of the parameters x + and x_. 

L-L(O) L-L(0) 




Figure 8. The top left-hand plot shows L(t) for the near extremal black hole, where there is no transition. 
The other plots show the case of first order behaviour discussed in the text, which exists far enough from 
extremality. The black hole is further from extremality as the plots are read from top left to top right, to 
bottom left to bottom right. 

Suppose at the energy E crit the same branch of geodesies has the least action. In that case 
there will be no source of non-analyticity at t crit , as the length of the geodesic (and therefore 
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the value of the correlator) will be a smooth function of the boundary point. It may then be 
that at some lower value of t one of the other branches will have the least action. This will 
generally result in non-analyticity similar to a second (or higher) order phase transition - 
the value of the correlator will be continuous past the transition, but some (time) derivative 
thereof will be discontinuous. 

On the other hand, behaviour analogous to a first order transition will result if one of 
the branches which develops at t CT n immediately dominates the correlation function. There 
is then no reason for the lengths of the geodesies on both branches to be degenerate, and 
in general one will discover a discontinuity of the correlation function, as opposed to its 
derivatives. 

For the values of the parameters we plot in figure 8, we see that when the transition 
exists it is always of first order, yielding a discontinuity in the correlator. By continuity it is 
clear that by tuning parameters we can make the two branches degenerate; thus as we vary 
the black hole parameters we obtain a line of first order transition terminating at a point of 
second order transition. 

4 Charged Correlation Functions 

We now discuss correlation functions of charged operators. As the analytic expressions are 
significantly more involved, we confine ourselves to discussing the qualitative features of the 
correlators, which are fairly similar to the neutral case discussed above. 

To approximate correlators of operators which carry R-charge, we need to discuss tra- 
jectories of electrically charged particles 13 . We will only consider radial motion here, the 
addition of angular momentum having much the same effect as in the uncharged case 14 . The 
Lagrangian describing such a particle with charge q and unit mass is 

C = \g^X» + q A»X» = \ (-/(r)i 2 + ^y) - ^ (4-1) 

where we have defined Q = cqQ. 

The worldline Hamiltonian 7i = (—f{r)t 2 + f~ 1 (r)r 2 )/2 is conserved, since the action is 

time translation invariant on the worldline. As we are interested in spacelike paths we can 

13 Such particles exhibit interesting phenomena such as Schwinger pair production and induced emis- 
sion See, e.g., and, for a recent discussion, PSj. The Schwinger effect in curved space is studied 
in EU. 

14 As in that case, for an arbitrarily small angular momentum, one can find trajectories probing the region 
beyond the Cauchy horizon, but they return to a copy of the boundary on the same side of the Penrose 
diagram. These are therefore irrelevant to the two-sided correlators discussed here, and more generally we 
expect such trajectories not to dominate any gauge theory correlator. 
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rescale A to set Ji = 1/2. The equations of motion are then 



t = 



f(r) 



1 



( 



r 2 ) ' 



( 



(4.2) 



Note that shifting A t by a constant (denoted in section 1 by $), which is just a gauge 
transformation, uniformly shifts the energy of all trajectories. This gauge freedom was fixed 
above, so that A t vanishes at the outer horizon. We choose here to absorb the constant part 
of A t , into the definition of E. 

The relations following from (j4.2j) are as follows. First, the (real part of) the boundary 
time, for symmetric trajectories, is 



where is the turning point described below. The length of the corresponding trajectory 
is given by 



where r max is the long distance bulk cutoff introduced in section 3. 

The equation for r(A) given in ()4.2|) cannot be interpreted as describing a particle with 
energy E 2 moving in a potential V(r), since now the potential itself depends on E. However, 
one can still describe the qualitative features of the trajectory, which now depend on both 
the parameters Q and E. 

The turning point, ve, of the trajectory is the largest root of the equation 



To visualise the situation, in figure 9 we plot V(r) = —f(r) and the (positive semi-definite) 
function (E — Q/r 2 ) 2 , the behaviour of which depends on the sign of Q/E. We therefore 
plot the function for both signs. 

In either case, there are three possibilities: 

• The graphs do not intersect, and there is no classical turning point. These classical 
trajectories end at the singularity. 

• The graphs intersect twice, in which case the turning point is that intersection with the 
largest value of r. As in the uncharged case, this will always be in the region between 
the inner and outer horizons. 




(4.3) 




(4.4) 




(4.5) 



21 



Figure 9. The left-hand plot has Q/E < 0. V(r) = -f(r) is in red and (E - Q/r 2 ) 2 in blue. Three 
trajectories with the same energy, E 2 , but different charges, are shown. The right-hand plot has Q/E > 0, 
and shows four trajectories with various values of the parameters. One has Q/E > without a radial 
turning point, one has Q/E = ri. and two have Q/E < r_, one with a radial turning point and one without. 

• The borderline case, in which there is a single intersection point, now not necessar- 
ily at precisely the top of the potential. Similar to the discussion above, looking at 
the behaviour of (J4.3|) and (J4.4j) near te reveals logarithmic divergences both in the 
boundary time and in the proper proper length of the geodesies. Therefore this value 
of E should be viewed as a limiting value, trajectories with larger value of E being 
irrelevant for our purposes. 

These features of r(A) are similar to those discussed previously. It is clear that only 
those trajectories with a radial turning point will return to a boundary, but it is not clear 
which boundary this will be. To determine which trajectories are relevant for our two-sided 
correlators, we must also consider the behaviour of t(X), which is more involved than in the 
uncharged case. Let us take E > so that t increases with affine parameter from the initial 
boundary 15 . Then, for Q < (the particle having a charge of opposite sign to the black 
hole), t just depends on the sign of the function f(r): for r > r + it is positive, whereas for 
r_ < r < r + it is negative. So, if such trajectories have a radial turning point, and they do 
for small enough Q, then they will connect the two boundaries of interest. In this case the 
coordinate t(X) is similar to the neutral geodesies. 

The situation for Q > (the particle having a charge of the same sign as the black hole), 

is more complicated. The behaviour depends on the magnitude of Q/E, and in each case is 

given by inspection of the right-hand side of the equation for t(X) in (|4.2J) : 

15 Trajectories with E < 0, but the same sign of Q/E, will just be mirror images of the trajectories we 
consider. 
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• For Q/E > r+, there is some r > r + for which i = 0, and the trajectory will start 
to move backward in coordinate time before reaching the outer horizon. It crosses 
the event horizon, upon which it moves forward in coordinate time again. If such 
trajectories have a radial turning point (and some certainly do), then they will connect 
the two boundaries of interest. 

• For Q/E = r^_, t — precisely at the outer horizon, where the particle will neces- 
sarily have a radial turning point. It will then return to the same boundary, so such 
trajectories are of no concern to us. 

• For r 2 _ < Q/E < r+, t does not vanish unless r_ < r < r + . However, it necessarily 
has a radial turning point before the radius at which t = 0, so will indeed connect the 
two boundaries of interest. 

• Finally, for Q/E < r 2 _, the particle may or may not have a radial turning point. If it 
does, then it will connect the boundaries of interest. Otherwise, the sign of t will flip 
at some r < r_, and the trajectory will always hit the singularity. 

The qualitative behaviour of the charged trajectories which do connect the two boundaries 
is thus similar to the neutral spacelike geodesies analysed above. They never cross the Cauchy 
horizon and, again, only reach a radius strictly larger than r_. 

5 Scanning Behind the Horizon 

5.1 Perturbing the AdS-RN Spacetime 

When the charged black hole is perturbed, linearised analysis suggests strong back-reaction 
near the inner horizon PJ El IS] • This is one instance of physics behind the horizon which 
would be nice to interpret in the gauge theory. We take here some preliminary steps towards 
this goal, postponing a more complete discussion for j26j . 

First, if one restricts attention to null perturbations which are either purely ingoing, 
or purely outgoing 16 (i.e., they are chiral in the sense of [HE]), then one can solve for the 
perturbed metric exactly EH] • One obtains an AdS version of the charged Vaidya so- 
lution [3^1 E2]- This solution contains one arbitrary function (the mass function) which 
depends on the profile of the perturbation. That profile can be varied, and can be arbitrarily 

16 For a massless field in AdS space, there are two types of perturbations, the non-normalisable one goes 
to a constant near the boundary. These are the type of perturbations we consider here. In the gauge theory 
this corresponds to turning on a marginal perturbation, with a specific time-dependent profile. 
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localised 17 . By investigating the two-sided correlators (and other observables) as a function 
of the perturbation, one obtains an efficient method of scanning (at least some of) the region 
behind the horizon. We will demonstrate here the effects on our observables of a simple wave 
profile, that of an infinitely localised delta function pulse. 

Secondly, in the asymptotically flat case, when both types of null perturbations are turned 
on 18 the spacetime is drastically changed near the inner horizon, resulting in the phenomenon 
of mass inflation j2H|. We expect that this phenomenon persists in the asymptotically AdS 
case jSZl I2E]- We further expect that the dependence of some gauge theory observables 
will be non-analytic as a function of the perturbation strength: taking that strength to be 
arbitrarily weak does not diminish its effect. 

Such non-analytic behaviour is not uncommon in field theories, and results from the 
existence of infra-red divergences. When re-summing perturbation theory, one discovers 
non-analytic dependence on coupling constants. Such behaviour would be slightly surprising 
for gauge theory on a compact space (a three-sphere in this case), at finite temperature and 
density, but it may be sensible in the infinite N limit. 

Here, we find a behaviour which is perhaps more sensible. As we saw above the spacelike 
geodesies accumulate at some finite distance from the Cauchy horizon. This is no longer 
necessarily true of the perturbed geometry. Nevertheless the geodesies we investigate seem 
to always be screened from the dramatic behaviour at the Cauchy horizon, as is demonstrated 
below. It is unclear if this is a property of the particular observables we are studying, or a 
more general property of the gauge theory. 

An Asymptotically AdS Charged Vaidya Solution 

The Vaidya solution [3^] describes an ingoing or outgoing null shell of matter incident on a 
vacuum black hole. Poisson and Israel generalised this to include charge in (22], their aim 
being to investigate the properties of the Cauchy horizon of the asymptotically flat Reissner- 
Nordstrom black hole. It is easy to further generalise these solutions to the asymptotically 
AdS case. In fact, though we will not do this here, one can derive J2.T3.I asymptotically AdS 
versions of the solutions relevant to mass inflation [2H] , describing both ingoing and outgoing 
flux. 

In terms of the lightcone coordinate u = t — r* defined in section 2, the outgoing solution, 
in the region r_ < r < r + , has the metric 

ds 2 = -f(u, r)du 2 - 2dwdr + r 2 dQ 2 3 (5.1) 

17 This arbitrariness is similar to gravitational plane waves which solve Einstein's equation for an arbitrary 
wave profile. 

18 These are still not generic perturbations, so cosmic censorship considerations do not apply. 
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where all the ^-dependence in the metric function is through that of the mass function M{u): 

M(u) Q 2 „ . x 

/Kr) = l--^ + ^ + r 2 . (5.2) 

The ^-dependence is set by the flux sent in from the left-hand boundary, which must have 
the null energy-momentum tensor, whose only non-zero component is 

ld u M(u) 



■ uu 



2 r 3 



(5.3) 



It is easy to check that this metric and energy-momentum tensor, together with the back- 
ground gauge field given in (j2.1|) solve the five dimensional Einstein-Maxwell field equations 
with negative cosmological constant. 

A Simple Example 

The simplest case, which we will consider here, is to take 

M(u) =m 1 + Q(u - u p ) Am. (5.4) 
The corresponding energy-momentum tensor 

r- - (5-5) 

corresponds to an infinitely thin shell of null matter, the flux leaving the left-hand boundary 
at t — t p . The spacetime splits into two regions: for u < u p the mass is mi and for u > u p , 
the mass is m2 = mi + Am. We can then apply the previous analysis to the two regions, 
each with its own potential, and match across the u = u v surface. For non-rotating neutral 
geodesies the potential in each region is simply V(r) = -/ mi;! (r), where / m (r) is the radial 
function ()5.2j) with specific value of the mass parameter m. 

If we take Am > 0, then f m2 < f mi for all values of r. The roots X-(m2) and x + {m2) 
of f m2 will thus be respectively less than i_(mi) and greater than x + (mi). The resulting 
geometry is shown schematically in figure 10. 

Now, to see the effect of the perturbation on the gauge theory, consider again the two- 
sided correlator in the geodesic approximation. The perturbation is then encoded in prop- 
erties of spacelike geodesies which intersect the surface u = u p . Some possible geodesies are 
shown schematically in figure 10, in which we denote the intersection point by ri,tj. The 
interaction with the flux will change the path of the geodesic beyond this intersection point. 
We are only interested in those geodesies which intersect the flux in the region beyond the 
outer horizon and it is clear that late-time correlators will be most suitable. 
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Inner apparent horizon 




Event horizon 



Figure 10. The geometry of the perturbed spacetime, showing the event horizon, the Cauchy horizon and 
the inner apparent horizon. The perturbation leaves the left-hand boundary at t = t p , and is shown as a 
solid red line. Spacelike geodesies are shown schematically in blue, the point at which they interact with the 
perturbation being shown as r/,t/. 

As we change the time t p at which the pulse leaves the boundary, the intersection point 
will change, and so will the length L of any geodesic. We should thus be able to see how 
the correlator changes as a function of t p , a clear signal of physics behind the horizon. We 
demonstrate here how the information about the perturbation can be detected in our simple 
example. 

The only subtlety 19 is that the energy of the geodesies is no longer conserved as they 
cross the null shell. In the first region, from the metric (|5.1J) . the conserved quantity is 

1 dC 

E = ~2dil = f ( r ^ + f (5-6) 

which is equal to the conserved energy E = — ||£ in the t, r coordinates we have been using 
for the unperturbed geometry. In the perturbed case, however, E is no longer conserved: if 
we define E = f(r,u)u + f then, from the geodesic equations, 

^ = ld u f(r, U y (5.7) 

19 As pointed out to us by Simon Ross. 
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which, for the mass function as in (|5.4|) . gives 



dE 



1 Am 



5{u — Mp) iidu. 



(5.8) 



2 r 2 



Whereas the "energy" E is discontinuous across the null shell, the momentum normal to 
the shell, which is just p± = it, is continuous at the interaction point jHHl- We can thus 
integrate the above equation, taking u = u(mi), the value of the momentum in the region 
with mass 7%. From the definition of E, and using the Lagrangian constraint C = 1 for 
spacelike geodesies, we have 



where E\ and E 2 denote the respective energies in the regions with mass m\ and rri2. Note 
that f(rj,u p ) < so that E 2 > E\. 

5.2 Detecting the Perturbation 

Qualitative Features 

To analyse the qualitative properties of the geodesies, we re-consider the particle mechanics 
problem of section 3. The particle moves from the right-hand asymptotic region to the 
intersection point at r = r/ in the effective potential V mi {r) = —f mi (r), and with energy E\. 
After the intersection point, the particle continues on its motion, but now in the effective 
potential V m2 {r) = — / ma (r), and with energy E 2 . The various possibilities are shown in 
figure 11, in which we fix the energy of the geodesic, and vary the mass, which in turn 
changes the potential. We can vary the intersection point by varying t p . There are then 
three possibilities (we will assume that the particle continues its motion entirely in V m2 after 
the intersection, i.e. that it does not intersect the pulse a second time): 

• If the intersection point satisfies rj > r^rr^), the turning point of the motion in 
V ms (r), then the particle continues its motion in the modified potential until reaching 
the new turning point r^^), and will thus still return to the opposite boundary. This 
is shown in the left-hand plot in figure 11. 

• The centre plot has 77 < r^^), so that after the interaction, it seems that the 
particle would be moving below V m2 (r). However, the new energy in the second region 




(5.9) 



so that 




(5.10) 
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is E 2 > Ei, and this shift in the energy across the null pulse is sufficient to cause 
the particle to jump back above the second potential (we have checked that the shift 
in energy is large enough for a range of black hole parameters). This case is thus 
qualitatively similar to the previous one. 

• The right-hand plot shows two ingoing particles, both with E\ > V mauX (rrii), so they 
would have fallen into the singularity in the unperturbed case. Here, however, they 
still have 77 > 775(7712), so will continue their motion until reaching this turning point. 
Contrary to the unperturbed case, such geodesies will thus return to the opposite 
boundary. 

The last case is interesting, since the trajectories now seem to cross the Cauchy horizon. 
However, as shown schematically in figure 10, by studying the trajectories more carefully one 
can see that this does not really happen. The trajectories cross the left-hand branch of the 
surface r_(mi), which would have been the Cauchy horizon in the unperturbed spacetime, 
but is not in the perturbed case. The geodesies never cross the surface r_(m 2 ), which is the 
left-hand branch of the Cauchy surface in this case. Once again, it seems that the correlators 
are geometrically protected from the catastrophic instability of the Cauchy horizon. 



V V V 




Figure 11. The particles with energy Ef are initially moving in the smaller potential, V mi (r). They interact 
with the perturbation at some r — ri < r+(mi), inside the outer horizon, at which point the energy jumps 
to E% > Ef, and they continue to move in the larger potential, V m2 (r). The various possibilities shown are 
discussed in the text. 

To demonstrate the scanning process of the region behind the event horizon, we will 
study the features of geodesies which fall into the first class discussed above, those with 
Ef < V max (mi) which intersect the flux at 77 > ^(7712). Such geodesies connect both 
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boundaries of the perturbed spacetime, and therefore dominate the correlation functions of 
the perturbed gauge theory. 



A Quantitative Plot of the Correlators 

It is difficult to solve analytically for the behaviour of the correlator in the perturbed space- 
time. Instead we demonstrate the process by numerically plotting the length of the geodesic 
as a function of t p , the time of the perturbation. 

Specifically , we plot here the fractional change in the length of the geodesic as a function 

L(Am = 0) K ' 

The method to generate the plot is not entirely obvious, so we outline it before presenting 
the result. 

We first fix the initial roots, x±(m,i), of the potential V mi (r), which fixes the initial mass 
mi and the charge (which remains unchanged). We fix Am, the strength of the perturbation, 
which fixes the final mass m 2 , and also the new roots a;±(m 2 ) of the potential V m2 (r). Finally, 
we fix the initial energy E\ of the geodesic in the unperturbed geometry, which fixes the 
boundary time of the correlator in that case. We take E\ ~ V max (mi) to give late time 
geodesies, which should be guaranteed to interact with the perturbation in the relevant 
region. 

In the unperturbed spacetime, with mass mi, the geodesic with energy E\ connects 
two specific points on opposite boundaries. The proper length of this approximates the 
correlator of operators inserted at these points. We need to compare this with the correlator 
in the perturbed spacetime, which connects the same two points. The correlators we compute 
before and after perturbing the spacetime should be the same; in the geodesic approximation, 
however, they are dominated by geodesies with different energies, say energy E\ in the first 
region of the perturbed geometry (and, as discussed above, with a shifted energy E\ in the 
second region). 

Let us denote the points on the left- and right-hand boundaries which the geodesies 
connect as U and t r . In the unperturbed spacetime, we have 

f°° dr 
t,(m0 + t r ( mi ) = 2E 1 / — __=_=_ (5.12) 

Jr El ( mi ) V mi {r)y/E( - V mi (r) 

which can be computed as in subsection 3.2. This must be equal to the same quantity in 
the perturbed spacetime: 

~ [°° dr ~ f ri dr 

ti(m 2 ) + t r (mi) = E 1 j = - E 1 



V mi (r)JE* - V mi (r) J W m ^ V mi (r)JE* - V mi (r) 
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"°° dr - C ri dr 
+ E 2 I + E 2 , (5.13) 



V m2 (r)^El - V m2 (r) J ^ 2 (m 2 ) V m2 (r)y/El - V m2 (r) 

as can be seen by inspection of figure 10, and where E 2 is related to E\ as in ()5.10|) . Again, 
we can compute this in terms of elliptic integrals as in subsection 3.2. In the perturbed case 
the result is a function of the intersection radius rj. Equating the two expressions for the 
boundary time ()5. 1215. 13|) in principle gives then a relation rj(Ei). In practice, however, we 
need to know the interaction radius to compute E 2 (through (j5.10|0 . which in turn we need 
to compute the interaction radius. 

We can avoid this vicious circle by running through values of rj (knowing that r + (mi) > 
rj > r^(m 2 )) for each E\, computing E 2 at each step (we take the smaller solution, with 
the negative sign in ()5.10|) . to ensure E\ < V mauX (rrii)) , then comparing the two expressions 
(I5.12I5.13|) for these specific values. If they are close enough, then we know that we have the 
correct value for rj. Numerical experiment for various values of the black hole parameters 
shows that the difference between the two expressions (j5.12l5.13j) for the boundary time is 
initially negative, for rj ~ r + (mi), and goes positive for rj ~ r_(m 2 ). So we are guaranteed 
a value of rj for which the two expressions are identical, and we simply take that value for 
which they are closest. 

We can use the resulting relation ri(Ei) to plot the proper length, L, and the time at 
which the spacetime is perturbed, t p , as functions of E\. The proper length of the geodesic 
in the perturbed spacetime is 

TlK , f°° dr n dr 

L(Am) 



^(mi) ^E\ - V mi (r) Jr^ mi ) ^ E 2 _ y mi{r) 

/ , + r ~ (5 - 14) 

•«» W \E 2 2 - V m2 (r) J r &2 (™>) J El - V m2 (r) 



and the time t p is 



t p = E 1 j — — — — — = == — E x 



T E X (mi) ( r ) Jr &1 (mi) V mi (r) y/Ef- V mi (r) 

r ri dr 
+ Ei -r.(rj) +n(oo). (5.15) 

Jr ^V mi {r)yjEl-V mi {r) 

Both expressions can, again, be evaluated in terms of elliptic integrals as in section 3. Both 
L and t p are functions of E\, utilising the relation n(Ei) found previously. 

To recap, we run through values of Ei, then run through values of rj at each step, 
computing E 2 , and evaluating (|5.12|) and ()5.13|) . We find a relation ti(E\) by taking those 
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values of 77 for which ()5.12|) and (j5.13j) are approximately equal, then evaluate L and t p 
for each {Ei,ri}. It is then a simple matter to plot the function 5L(t p ) as given in (|5.11j) . 
The result is shown in figure 12. This schematic plot is for some specific values of the 
parameters (x_ ~ 0.25, x + ~ 2.5, which gives mi = 10 and Q = 1.5, Am = 0.5 and 
E 2 ~ 0.999V max (m 1 )), but it seems to be representative of the general behaviour. 

5L 




tp 

Figure 12. Here we show the fractional change in the length of the geodesic as a function of t p , for some 
choice of parameters as explained in the text. The lack of smoothness for small t p is due to numerical effects, 
and is not expected to be physical. 

This plot demonstrates that our correlators are sensitive to a perturbation of the geometry 
which is localised beyond the outer horizon. We postpone further discussion for future work, 
and conclude by making a few comments. 

First, we should note that figure 12 cannot be used for the region V max (mi) < E\ < 
V max (m2) of energies, since in this case the integrals have to be evaluated differently. Similar 
methods should work to generate a plot of SL(t p ) for that range of energies however. 

Moreover, we can only probe part of the region of the spacetime beyond the event horizon. 
There are no real solutions for 77 beyond a certain range of energies, and this range does not 
cover the entire region behind the event horizon. 

Finally, it is intriguing that the plot of L(t p ) in figure 12 is not single- valued. This means 
that there are two geodesies, with different energies E\, which connect the fixed boundary 
points in the perturbed spacetime. Similar to the case of the unperturbed spacetime, this 
is reflected in some non-analytic behaviour of the correlators in the perturbed spacetime. It 
would be nice to get a better handle, through analytic calculations, on these issues. 
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Appendix: The t(E) Integral 

Consider the integral ()3.5|) for non-rotating geodesies. Then V(r) = — /(r) and, dropping 
the residue coming from integrating over the horizon, we have 

E f°° x 5/2 

-t = — dx = (5.16) 

2 Jx + (x + x )(x — X-)(x — x + )y/(x + Xq){x — xJ){x — x + ) 

Defining a new variable t through 



with a a constant to be determined, gives 

77. pa / \ -1/2 / ~ \ 5/2 

^ A I jj. / 1 C 0^2 \ X - .2 \ I ft C 0,2\ ( -i C -j-2\(i H-.2 



-t = -A dt 1 - -f] 1 - —t z 1 - -t z 1 - — t z 1 - -^t 

2 / n V a I V ax+ J \\ a / V a J \ a 



(5.18) 



where we have set 



A = 
a 



2x 5/2 



THJ- — "7 ~ m/2 ((^+ + x o)( x + ~ — x +)) 

L ' 2 {X+ + Xq) l/Z 



X— "T Xq X— "T Xq X— — X— X— — , , 

c = - — — — , c = - — ■ , c_ = , c+ = . (5.19) 

Now observe that 

i - C ±A (i - ^; 2 ) (i - ^r 1 = co (i - ^"V (i - c -e\\c + (i - C -±A 1 

a J \ a J \ a J J \ a J \ a J \ a 

(5.20) 



where 

.2 „2 

c 



Cq a c_ a c 4 



(c - c_)(c - c+)' (c^ - Cq)(c^ - c + Y r (c+ - Co)(c + - c_)' 

(5.21) 
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If we take 



2 x- 1 

a = c , fc = , 

x + a 



a 



a — 



a 



Ot, 



(5.22) 



so that 

c 



ac 



(ckq — a?.)(ao — a +) ' 
then the integral becomes 



(a?. — «o)(«^ — « 2 1 



(a 2 — al)(a% — a 2 ) ' 
(5.23) 



to = —A I dt 



(1 - k 2 t 2 y 



c 



+ 



+ 



c+ 



1-t 2 V(l-«gt 2 ) (l-a 2 * 2 ) (l-a^t 2 );' 



(5.24) 



We now use 
'1 -k 2 t 2 1 



(a 2 - k 2 



1-t 2 (1-aH 2 ) 
and, noting that 



+ 



k 2 



(1 - aH 2 ) ^(l-t 2 )(l-k 2 t 2 ) a 2 y/(l-t 2 )(l-kH 2 ) 

(5.25) 

(5.26) 



^ + — + ^ = 

9 1 9 1 9 w 5 

ai a+ 



we have 



where 



E 



t n = —A / dt 



(1 - k 2 t 



2+2\2 



2 Jo v /(l-t2)(l_ fc 2 t 2 ) 

«o(«o - fc2 ) 



(Jo + J_ + J 4 



(qq — a 2 )(qq — a+) (1 — apt 2 )' 

q; 2 (q; 2_ A .2N j 



(a 2 - a 2 ) (a 2 - a 2 ,) (l-a 2 .* 2 )' 



ai(a; 2 



K-a^K-ai) (l-c^t 2 )' 
Each term can now be evaluated [2E]> to give 



(5.27) 

(5.28) 
(5.29) 
(5.30) 



+ 



-t = —A 
° 2 



(a 2 - k 2 f 



2 2 2 



■F(0,fc) + 



(a 2 - k 2 f 



-11(0, a 2 , A;) + 



«o(«q — a 2 _)(o;o — a+) 
(a 2 - A; 2 ) 3 



n(0,ao,fc) 



a 2 (a 2 — aQ)(a 2 — a+) vr '~ v ' a+(a+ — a 2 ,)^ 2 — a 2 ) VY ''""+' 7 J ^ 
where = sin -1 a. Substituting for the various constants gives the result ()3.7j) in the text. 
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